Graphene quantum dots in perpendicular magnetic fields 
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We report transport experiments on graphene quantum dots. We focus on excited state spectra 
in the near vicinity of the charge neutrality point and signatures of the electron-hole crossover 
as a function of a perpendicular magnetic field. Coulomb blockade resonances of a 50 nm wide 
and 80 nm long dot are visible at all gate voltages across the transport gap ranging from hole to 
electron transport. The magnetic field dependence of more than 40 states as a function of the back 
gate voltage can be interpreted in terms of the unique evolution of the diamagnetic spectrum of 
a graphene dot including the formation of the E — Landau level, situated in the center of the 
transport gap, and marking the electron- hole crossover. 

PACS numbers: 73.22.-f, 72.80.Rj, 73.21. La, 75.70.Ak 



INTRODUCTION 



Graphene nanostructures 0, 0, 0, |, 0, 0, 0, 0, SEE] 
attract increasing attention mainly due to potential ap- 
plications in high mobility electronics Il2lll3j and solid 
state quantum information processing [14j |. In particu- 
lar, low nuclear spin concentrations expected in graphene 
promise long spin lifetimes [3, fl5. \uL Q3] and make 
graphene quantum dots (QD s) , H, 0, 0, 0] interest- 
ing for spin-qubit operations I Li]. Moreover, graphene 
nanostructures may allow to investigate phenomena re- 
lated to massless Dirac Fermions in confined dimen- 
sions [H, [ljj 0, Intensive research has 
been triggered by the unique electronic properties of 
graphene [23[ including the gapless linear dispersion, and 
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the unique Landau level (LL) spectrum 
search for signatures of graphene-specific properties in 
quantum dots is of interest in order to understand the 
addition spectra, the spin states and dynamics of con- 
fined graphene quasi-particles. Recent advances in fab- 
ricating width-modulated graphene nanoribbons helped 
to overcome intrinsic difficulties in (i) creating tunneling 
barriers and (ii) confining electrons in graphene, where 
transport is dominated by Klein tunneling-related phe- 
nomena [H, [2?| . Graphene quantum dots have been fab- 
ricated and Coulomb blockade [4 [5[, quantum confine- 
ment @, 0] and charge detection @ have been observed. 

Here, we show tunneling spectroscopy (i.e. transport) 
measurements on a fully tunable graphene quantum dot. 
We present the evolution of a large number of Coulomb 
resonances near the charge neutrality point in a magnetic 
field from the low-field regime to the regime of Landau 
level formation. In particular, we investigate the quan- 
tum dot spectrum in the vicinity of the charge neutral- 
ity point as a function of a perpendicular magnetic field 
(B-field). Near the electron- hole crossover we observe 
a rich excited state spectrum with inelastic cotunneling 
lines and unconventional features inside the Coulomb di- 
amonds. Some of these features, which are very pro- 
nounced at B = 0, quickly disappear in a perpendicu- 
lar magnetic field. Moreover, we find indications of the 
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FIG. 1: Low bias (V& = 100 /iV) back gate characteristic at 
V pg = 2 V. The resolved transport gap (see highlighted region) 
separates between hole and electron transport. This trace is a 
back gate cross-section of Fig. 1(e) in Ref. 0, which includes 
the region [B] marked therein. The dashed line marks the po- 
sition of the electron- hole crossover identified in Ref. [H . The 
inset shows a scanning force microscope image of an etched 
graphene quantum dot device with source (S) and drain (D) 
leads and a plunger gate (PG) for electrostatic tunability. The 
scale bar corresponds to 100 nm. 



formation of the lowest (E = 0) Landau level at high 
B-fields marking the electron-hole crossover in graphene 
devices. 

This paper is organized as follows. In section 2 we 
briefly review the graphene quantum dot device fabrica- 
tion. Transport measurements at zero magnetic field are 
discussed in section 3.1. Landau level formation and the 
magnetic field dependence of inelastic cotunneling pro- 
cesses are addressed in sections 3.2 and 3.3, respectively. 



DEVICE FABRICATION 

The state-of-the-art fabrication process of graphene 
nano devices, which has been mostly developed in Manch- 
ester liil [11, and in New York [251, is based on the 



2 




4.75 4.8 4.85 4.9 4.95 



Plunger gate V pg (V) 



FIG. 2: (color online) (a) Coulomb blockade resonances of the 
graphene quantum dot as function of plunger gate voltage 



Vpg at Vba 



-0.9 V. The conductance has been measured 



at constant bias of Vb = 100 /iV. (b) Coulomb diamonds, 
i.e., finite-bias measurement of the quantum dots differential 
conductance G q d — dl/dVb as function of Vb and V pg at B = 
0. Note the rich excited state spectra and in particular the 
different vertical distances between the diamond edge and the 
diagonal structures inside the Coulomb-blockaded region for 
positive Vb. 



mechanical exfoliation of (natural) graphite by adhesive 
tapes (28|. The substrate material consists of highly 
doped silicon (Si ++ ) bulk material covered with approx- 
imately 300 nm of silicon oxide (S1O2), where thickness 
and roughness of the Si02 top layer is crucial for the iden- 
tification and further processing of single-layer graphene 
samples. Standard photolithography followed by met- 
alization (chrome/gold) and lift-off is used to pattern 
arrays of reference alignment markers on the substrate 
that are later used to re-identify locations (of individual 
graphene flakes) on the chip and to align further process- 
ing patterns. 

The graphene flake has to be structured to submicron 
dimensions in order to fulfill the device design require- 
ments (see e.g., figure in the abstract). We use a tech- 
nique based on resist spin coating, electron beam lithog- 
raphy (EBL), development and subsequent etching of the 
unprotected graphene. Two successive processing steps 
have been used including small modifications to decrease 
the minimum feature size of the graphene nanodevice. 
First we used a resist [Polymethylmethacrylat (PMMA)] 
with a thickness of 45 nm and short etching time to define 
the delicate structures. In the second step a 100 nm thick 
resist is used for the coarse structuring of the flake also 
providing a broader process parameter window for the 



etching step. It has been found that short (5 and 15 sec- 
onds) mainly physical reactive ion etching (RIE) based on 
argon and oxygen (80/20) provides good results without 
influencing the overall quality of the graphene flake (29[. 
After etching and removing the residual EBL resist, the 
graphene nanostructures are contacted by an additional 
EBL step, followed by metalization and lift-off. Here we 
evaporated 2 nm chrome (Cr) and 40 nm gold (Au) for 
contacting the graphene quantum dot device. 

A scanning force microscope image of the device stud- 
ied here is shown as an inset in Fig. 1. The 50 nm wide 
and 80 nm long graphene quantum dot is connected to 
source (S) and drain (D) via two graphene constrictions 
with a width of 25 nm, both acting as tunneling barriers. 
The dot and the leads can be further tuned by the highly 
doped silicon substrate used as a back gate (BG) and an 
in-plane graphene plunger gate (PG). 



MEASUREMENTS 

All measurements have been performed in a dilution 
refrigerator at a base temperature of T « 90 mK. We 
have measured the two-terminal conductance through the 
graphene quantum dot device by applying a symmetric 
DC bias voltage Vb while measuring the current through 
the quantum dot device with a noise level below 10 f A. 
For differential conductance measurements a small AC 
bias, Vb,ac — 40 [iSf has been superimposed on Vb and the 
differential conductance has been measured with lock-in 
techniques at a frequency of 18.9 Hz. 

In Fig. 1 we show the differential conductance as a 
function of back gate voltage at low bias (V& = 100 /iV) 
highlighting the strong suppression of the conductance 
around the charge neutrality poi nt ( < Vb g < 7.5 V) 
due to the so-called transport gap 30, Hl[ . Here we tune 



transport from the deep hole to the deep electron regime, 
as illustrated in Fig. 1. The large number of resonances 
within the gap region may be due to both, (i) resonances 
in the graphene constrictions acting as tunneling barri- 
ers [4[ (and thus being main responsible for the large 
extension of this transport gap) and (ii) Coulomb reso- 
nances of the quantum dot itself. 



Coulomb blockade measurements at B = 

By fixing the back gate voltage at a value close to 
the transport gap (e.g., Vb g = —0.9 V) and sweeping 
the lateral plunger gate voltage (V pg ) in a narrow range 
close to the charge neutrality point (see below) Coulomb 
blockade resonances of the graphene quantum dot can be 
well resolved, as shown in Fig. 2(a). The strong ampli- 
tude modulation of the conductance peaks [see e.g., on 
the very left of Fig. 2(a)] is mainly due to transparency 
modulations of the constrictions [4( , which can dominate 
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FIG. 3: Evolution of Coulomb blockade resonances in perpendicular magnetic field. The conductance is plotted as function of 
back gate voltage Vbg and magnetic field B (Vb = 500 /jlV). The effective plunger gate voltage V p * g , calculated via the relative 
lever arm, is given on the top of the graph. This data have been measured in the regime marked by the solid line [B] in Fig. 1(e) 
of Ref. [|. 



the width of the Coulomb resonances and significantly 
elevate the conductance. In Fig. 2(b) we show corre- 
sponding Coulomb diamond measurements [see black bar 
in Fig. 2(a)], that are measurements of the differential 
conductance (G = dl/dVb) as a function of bias voltage 
Vb and plunger gate voltage V pg . Note, that this mea- 
surement was recorded at a higher base temperature of 
T = 370 mK. The reason for this slightly higher tem- 
perature is the good visibility of the diagonal lines inside 
the diamond. The overall structure of the diamond is 
not changed significantly by increasing the temperature. 
The logarithmic plot of the Coulomb diamonds measured 
in the close vicinity of the charge neutrality point dis- 
plays a rich excited state spectrum. We observe lines 
of increased conductance outside the diamonds running 
parallel to the diamond edges, which are well aligned 
to inelastic cotunneling onsets visible as faint horizon- 
tal (constant bias) structures inside the diamond-shaped 
regions (see arrows). At the diamond boundary, the hor- 
izontal lines seamlessly join some of the most prominent 
diagonal lines in the non-blockaded region, allowing to 
extract characteristic excitation energies in the range of 
2 to 3 meV. 

In addition, we observe for positive bias Vb > V un- 
conventional features inside the Coulomb-blockaded re- 
gions (see dashed lines) consisting of diagonal lines run- 
ning parallel to the diamond edge. The fact that they 
have the same slope as the diamond edges suggests that 
they are connected to the alignment of an energy level 
with source (negative slope) or drain (positive slope). 
The vertical distance between the diagonal lines and the 
diamond edge is identical for lines with positive and neg- 
ative slope (see dashed line inside the diamond). It has 
been shown by Schleser et al. Q that these features 
are related to sequential tunneling through the quantum 
dot, occurring after it has been excited by an inelas- 
tic cotunneling process. In simple terms, these features 
can be attributed to a transport configuration where the 
tunneling-out rate from the excited state is significantly 
larger than the relaxation rate to the ground state, which 
might be directly related to strongly non-linear energy 



dependences of the tunneling barriers. Finally, we also 
observe regions of negative differential conductance (see 
red areas), which are not yet fully understood. 



Coulomb resonances as a function of a perpendicular 
magnetic field 

In Fig. 3 we show more than 40 Coulomb resonances 
as a function of a B-field perpendicular to the graphene 
plane. The measurement has been taken in the back gate 
voltage range from Vb g = -3.5 to 6 V (also highlighted by 
the vertical line [B] in Fig. 1 and Fig. 1(e) in Ref. [§]). 
Thus we certainly tune from hole to electron transport 
(see arrows in Fig. 1). The evolution of Coulomb reso- 
nances in (perpendicular) magnetic field shows signatures 
of the electron-hole crossover of the quantum dot states 
(see below and Ref. [§]). There is a common trend of 
resonances at lower back gate voltages (see, e.g., reso- 
nances at Vbg = V) to bend for increasing B-field to- 
wards higher energies (higher Vb g ). In contrast we find 
for higher back gate voltage the opposite trend (see, e.g., 
resonances and arrow at Vb g ~ 5.8 V), where resonances 
tend to shift to lower energies for increasing magnetic 
field. This overall pattern is disturbed by additional fea- 
tures such as localized states, regions of multi-dot behav- 
ior and strong amplitude modulations due to constriction 
resonances. For example, we observe a strongly elevated 
background below Vb g < -1 V, which can be attributed to 
the increase of the transparency of the tunneling barriers 
in agreement with the measurement presented in Fig. 1. 
Additionally, we observe between Vb g = -1 and V a 
weakly coupled state crossing the Coulomb resonances 
at finite magnetic fields. We interpret the weak mag- 
netic field dependence and low visibility in transport as 
a manifestation of a strongly localized state. Moreover, 
we see several level crossings and splittings in the region 
above Vb g = 4 V, which might be due to the presence of 
an additional quantum dot spontaneously forming in one 
of the two constrictions, and coupling to the gates with 
different lever arms 
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FIG. 4: (color online) (a) Evolution of inelastic cotunneling onsets inside the Coulomb-blockaded region in perpendicular 
magnetic field. The conductance is plotted as function of bias Vb and magnetic field B = to 2 T (stepped in units of 0.1 T) 
for constant plunger gate voltage V£ g = 4.66 V [see vertical line in panel (b)]. For clarity the individual traces have been 
offset by 0.2 xl0 _3 e 2 /h. The peaks related to inelastic cotunneling marked by arrows quickly vanish as function of B. (b,c) 
Coulomb diamond measurements at B = 6.4 T (b) and B = 13 T (c), respectively. These measurements are taken at constant 
back gate voltage Vb g -0.9 V and highlight both (i) single dot signatures at high magnetic fields and (ii) the strong suppression 
of excited state-related signatures at higher B-fields. 



Finally, in close analogy to Ref. [8[ where the same 
sample was investigated in a different parameter regime 
we observe "kink" features with increasing B-field. It 
has been found that with increasing magnetic field the 
levels feature a kink signifying filling factor v — 2 in 
the dot (see dashed lines and arrows in Fig. 3) before 
converging toward the E = Landau level, in ana logy 
to what has been observed in GaAs quantum dots [341 ] . 
In particular, the population of the lowest Landau level 
leads to a region around E = @, [32[ with a remark- 
ably ordered addition spectrum, which emerges for large 
magnetic fields (see also the analytical calculations for 
a circular dot 0, [lo[ ) . This feature is a unique conse- 
quence of the interplay of the unique Landau level spec- 
trum in graphene |24l [25[ and the carrier confinement 
due to the finite size of the system. For more details, in- 
cluding numerical simulations of graphene quantum dots 
and a quantitative discussion of the peak-to-peak spacing 
evolution in B-field we refer to Refs. [8l l32j . 



Coulomb diamonds at finite magnetic field 



In Fig. 4a we show the evolution of inelastic cotunnel- 
ing onsets inside the Coulomb-blockaded region as a func- 
tion of perpendicular magnetic field. The conductance is 
plotted as a function of bias voltage Vb and magnetic field 
at constant plunger gate voltage V pg = 4.66 V [see verti- 
cal line in Fig. 4(b)] and the B-field has been stepped by 
0.1 T from B = to 2 T. While at B = T we observe sev- 
eral peaks and steps (black and white arrows), at B = 2 T 
the trace is smoother with remaining steps at 1.4 mV, - 
1.3 mV and 4.5 mV (white arrows). The different B-field 
dependence of the features might be attributed to two 
origins. The fine structures disappearing at B = 0.15 T 
(black arrows) may arise from energy dependent fluctu- 
ations of the barrier transmission, while the steps (white 
arrows) are attributed to the onset of inelastic cotunnel- 
ing channels. This is supported by Figs. 4(b), (c), where 
we show Coulomb diamond measurements at B = 6.4 T 
and B = 13 T. These measurements are taken at the same 
constant back gate voltage V bg = -0.9 V as in Fig. 2(b) 
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and are plotted with the same color scale. Thus, we can 
directly compare Figs. 2(b), 4(b) and 4(c) simplified by 
the diamond labels I and II [see Figs. 4(b), (c); the dia- 
mond shown in Fig. 2(b) corresponds to diamond II]. The 
horizontal (constant bias) lines in the Coulomb blockaded 
region present at B = T are absent at higher B-fields. 
However, at B = 6.4 T cotunneling onsets and the un- 
conventional feature consisting of diagonal lines inside 
the diamond [see above and dashed lines in Fig. 4(b)] 
can still be observed, but in a different diamond than at 
B = T (I instead of II). At B = 13 T all excited-state 
related signatures inside the diamond are suppressed. In 
addition, also the negative differential conductance out- 
side the diamond is significantly reduced and might also 
be related to a decrease in the energy dependent fluctu- 
ation of the barrier transmission. However, from these 
measurements we can conclude that the single-dot be- 
havior of this etched graphene quantum dot device stays 
well preserved also for high magnetic fields. 



CONCLUSION 

We have performed detailed studies of transport 
through a graphene quantum dot in the vicinity of the 
charge neutrality point. The evolution of Coulomb res- 
onances in magnetic field showed signatures of Landau 
level formation. Indications for the crossing of filling fac- 
tor v = 2 are obtained by the observation of "kinks" in 
spectral lines before bending towards the charge neutral- 
ity point. Coulomb blockade diamonds at B = T near 
the electron-hole crossover show many excited states and 
cotunneling onsets accompanied by diagonal lines in the 
diamond attributed to cotunneling mediated sequential 
tunneling through an excited state. Many small scale 
lines inside the diamond disappear quickly (B = 1.5 T) 
in a magnetic field while more step like features (cotun- 
neling onsets) are still visible at B = 6.4 T. The presented 
measurements open the way for an in-depth exploration 
of the few charge carrier regime including addition spec- 
tra in graphene quantum dots. 
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